Color neutrality effects in the phase diagram of the PNJL model 



D. Gomez Dumm, 1 - 2 B D. B. Blaschke, 3 ' 4 '[j] A.G. Grunfeld* , 2 - 5 % and N. N. Scoccola 2 ' 5 ' 6 <0 

1 IFLP, CONICET — Dpto. de Fisica, Universidad Nacional de La Plata, C.C. 67, 1900 La Plata, Argentina 
e CONICET, Rivadavia 1917, 1033 Buenos Aires, Argentina 
3 Institute for Theoretical Physics, University of Wroclaw, Max Born place 9, 50204 Wroclaw, Poland 
J < Bogoliubov Laboratory of Theoretical Physics, JINR Dubna, Joliot- Curie Street 6, 14-1980 Dubna, Russia 
J Physics Department, Comision Nacional de Energia Atomica, 

Av. Libertador 8250, 1429 Buenos Aires, Argentina 
6 Universidad Favaloro, Solis 453, 1078 Buenos Aires, Argentina 
(Dated: December 22, 2008) 

The phase diagram of a two-flavor Polyakov loop Nambu— Jona-Lasinio model is analyzed impos- 
ing the constraint of color charge neutrality. The main effect of this constraint is a coexistence of 
the chiral symmetry breaking (xSB) and two-flavor superconducting phases. Additional effects are 
a shrinking of the x^B domain in the T — fx plane and a shift of the end point to lower temper- 
atures, but their quantitative importance is shadowed by the intrinsic uncertainties of the model. 
' The effects can be understood in view of the presence of a nonvanishing color chemical potential Us , 

which is introduced to compensate the color charge density p$ induced by a background color gauge 
^"*^ ■ mean field 03. At low temperatures and large chemical potentials the model exhibits a quarkyonic 

phase, which gets additional support from the diquark condensation. 
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CN . PACS numbers: 12.29.Fe, 21.65. Qr, 25.75.Nq 

O I. INTRODUCTION 

v : 
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i ' Constraints on the phase diagram of QCD under extreme conditions of high excitation (temperature) and com- 
pression (density) are of vital interest for large-scale experimental programmes with ultrarelativistic heavy-ion beams, 
' searching for signatures of the QCD phase transitions: chiral symmetry restoration and deconfinement. 

Under conditions of finite temperature T and small chemical potentials fi, which are probed in heavy ion collisions 
at SPS, RHIC and in the near future at LHC, Lattice QCD simulations have provided insight into the phase structure 
since the temperature-dependent change of the order parameters of the above phases, namely the chiral condensate 
(qq) and the traced Polyakov loop <&, have been determined [H, 0j & 0] • The peak positions of the susceptibilities 
derived from these calculations for both transitions are remarkably coincident at the same pseudocritical temperature 
T c . This might be understood in terms of general effective theory [f|. It is also remarkable that these results can be 
nicely described with an effective model of the NJL type, extended with the inclusion of a coupling to the Polyakov 
loop [1, 0, H, @ with a potential U($,T) that can be extracted from Lattice QCD simulations of the pressure in the 
pure gauge theory [13, [H| • 

. £^ \ The challenge for experiments as well as for theory is to extend this knowledge into the domain of finite baryon 
densities, where precursors of color superconductivity (pseudogap phase [Hj]) or even color superconducting quark 
matter phases themselves can occur (la . Il4j and new constraints from observations of compact star properties may 
apply [ll El- 
Most of the effective model analyses addressing the chiral phase transition in the QCD phase diagram have been 
performed with NJL type models which, however, lack quark confinement and lead to an onset of finite quark densities 
already at unphysically low temperatures of about 50 MeV [I7j . far below the deconfinement temperature of about 
200 MeV 0,11 • 

The straightforward step to implement the effects of confining forces and to suppress unphysical quark degrees 
of freedom is to study the phase diagram of quar k matter in the Polyakov loop NJL (PNJL) model. This has 
been carried out, e.g., in Refs. [ll|, [H, 0j| I20I l2ll. [22j]. In these works, however, the important question of color 
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neutrality of quark matter has not been considered. The inclusion of neutrality constraints in the PNJL model 
has been taken up recently in Ref. (23| . In that paper, however, only the high-density region for a quark chemical 
potential fi — 500 MeV has been considered, and light quark masses have been neglected. It has been found that the 
Polyakov-loop formulation together with neutrality constraints can give rise to an increase of the critical temperature 
for the color superconductivity transition, in deviation from the well-known BCS relationship. An investigation of the 
interrelation between chiral symmetry breaking (xSB) and color superconductivity, as well as a study of the effect of 
color neutrality without enforcing electric neutrality, is still missing. The aim of the present paper is to address these 
subjects, analyzing the QCD phase diagram in the framework of a generalization of the PNJL model that fulfills the 
constraint of color neutrality. 

The article is organized as follows. In Section II we introduce the formalism for enforcing color neutrality in the 
PNJL model by color chemical potentials. In Section III we discuss the numerical results for the phase diagram, 
first in the non-superconducting case, then including color superconductivity. In the final section we present our 
conclusions and discuss possible consequences for the phenomenology of next heavy-ion collision experiments. 

II. FORMALISM 

The Euclidean action for the two-flavor PNJL at temperature T is given by 

r r _ G - 

dr d 3 x (- i^D^ + m)ip - — [(ipipf + (^Ts^) 2 ] 

- y (V ; c«75r 2 A a V)(^c«75r 2 A a V) t + W($,T)j , (1) 

where ip is the Nf = 2 fermion doublet tp = [u d) T , rh = diag(m„, stands for the current quark mass matrix, and 
f3 = 1/T. For simplicity we consider the isospin symmetry limit, in which m u = = rfiQ. Charge conjugated fields 
in Eq. (TTJ) are defined by ipc — 7274 while f and A a , with a = 2,5,7, stand for Pauli and Gell-Mann matrices 
acting on flavor and color spaces, respectively. This Euclidean action leads to local chiral invariant current-current 
interactions in the quark-antiquark and quark-quark channels. The latter is expected to be responsible for the presence 
of a color superconducting phase in the region of low temperatures and moderate chemical potentials. 

The coupling of fermions to the Polyakov loop has been implemented in Eq. ([T]) through the covariant derivative in 
the fermion kinetic term j^D^, where = d^ — iA^. Here the Euclidean operator 7^9^ is defined as 74^+7- V, with 
74 = *7o- As usual, we assume that the quarks move in a background gauge field A = g S^o G%\ a /2, where G% are the 
SU(3) color gauge fields. Then, at the mean field level, the traced Polyakov loop is given by $ = |Tr exp(i,00), with 
<t> = iA{) = constant. The traced Polyakov loop can be taken as an order parameter for the confinement transition. 
In the pure glue theory it can be associated with the spontaneous breaking of the global Z 3 center symmetry of color 
SU(3), <& = corresponding to the symmetric, confined phase [§4j]. Here we will work in the so-called Polyakov gauge, 
in which the matrix is given a diagonal representation — 03 A3 which leaves only two independent variables, 

03 and 08- Finally, the action ([1]) also includes an effective potential U that accounts for gauge field self- interactions. 

In order to obtain the grand canonical thermodynamical potential of the PNJL model at temperature T and chemical 
potential /i = /is/3 we start from Eq.{T]), performing a standard bosonization of the theory. This can be done by 
introducing fields (a, 7?) and A a corresponding to the quark-antiquark and diquark channels respectively. Then we 
consider the mean field approximation, keeping the nonzero vacuum expectation values a and <f> and dropping the 
respective fluctuations. Concerning the superconducting vacuum, we adopt here the usual ansatz in which one has 
A5 = A7 = 0, A2 = A. In this way, using the standard Matsubara formalism one obtains 

T °° f d 3 n a 2 A 2 

Omfa = - g E / 7^3 ln det + 2~G + 2H + m ' T) ' (2) 

where the inverse propagator S^ 1 is a 48 x 48 matrix in Dirac, flavor, color and Nambu-Gorkov spaces. The dependence 
of S~ x with T and [x (fj, = /^s/3) is obtained from the T = /i = four-momentum integrals in the fermion determinant 
by replacing P4 — > uj n — iji, where oj n = (2n + 1)ttT are the usual fermionic Matsubara frequencies, and fi stands for 
a complex "chemical potential" matrix in color space, 

jii = diag(/x r , fj, g ,^b) + « (<fe A 3 + 0s A 8 ) , (3) 

with 



diag(/i r ,/i s ,^ h ) = fit + /13A3 + /i 8 A 8 



(4) 
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In general, the inclusion of color chemical potentials ^3 and /ig is necessary to ensure color neutrality. Concerning 
the mean field effective potential £/($,T), which accounts for the Polyakov loop dynamics, we use here a form that 
appears to be consistent with group theory constraints as well as lattice results (from which one can estimate the 
temperature dependence). Following Ref. [ll| we take the logarithmic form of the Polyakov-loop potential, motivated 
by the SU(3) Haar measure 

= -ia(T)<r$ + b(T) ln[l-6$*$ + 4($ 3 + $* 3 )-3(<r$) 2 ] , (5) 

with the corresponding definitions of a(T) and b(T) from 

It is worth to notice that in presence of the Polyakov loop the thermodynamical potential could be in general a 
complex quantity. As discussed in Ref. jl9j, in order to properly define the mean field approximation one should 
require that the mean field configuration provides the maximal contribution to the partition function. Then, the 
mean field values (order parameters) a, A, 3 and <j)8 should fulfill the coupled set of "gap equations" 

<9Re [^mfa] „ 

- 0. (6) 



(da, dA, a03,50 8 ) 

If one also imposes color charge neutrality conditions, one has two additional equations that fix the color chemical 
potentials /X3 and /is- The vanishing of color charges implies 

<9Rc [Omfa] = Q ^ 
(<9/u 3 ,<9/i 8 ) 

Finally, in order to obtain a physically meaningful solution, it has to be verified that the resulting field configuration 
leads to a real-valued thermodynamical potential. 

At vanishing chemical potential, owing to the charge conjugation properties of the QCD Lagrangian, the mean field 
traced Polyakov loop $ is expected to be a real quantity. Since 03 and (f>s have to be real- valued (l9j . this condition 
implies 0§ = 0, and it is easy to see that the thermodynamical potential turns out to be real, as desired. In general, 
this need not to be the case at finite [i [25ll26l [27j. As in Refs.fTj], [HI [23[ we will assume that the potential U is such 
that the condition 0g = is well satisfied for the range of values of /i and T investigated here. The traced Polyakov 
loop is then given by <f> = <f>* = [1 + 2 cos(/303)]/3. In addition, it can be seen that if 0§ = then Eqs. ([6]) and ([7]) 
lead to /i3 = 0, leaving /i 8 as the only potentially nonvanishing color chemical potential. This also implies that the 
condition = is trivially satisfied. Notice that the fact that /13 = 0, together with the ansatz chosen for A a , allow 
for a residual r — g color symmetry, leading to fi r = 

The Matsubara sums can be now explicitly evaluated. One obtains 
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^mfa(T) = -2 / ^^[2T ln(l + e-^) + E s ] + ^ + ^ + W(*,T) , (*) 



where the quasiparticle energies Ej are given by 



-El, 2 


— EpTW 




-^3,4 




-A 2 


-^5,6 


= V ( £ p T Ms) 2 " 


h A 2 



(9) 

z0 3 (10) 

Z0 3 ■ (11) 



Here we have defined e p = \/\p\ 2 + M 2 , being M — ttiq + a the quark constituent mass. 

The momentum integral is of course divergent. Following the usual NJL regularization prescription, we cut the 
integrals of the zero-point energies Ej at a given cutoff A, while the logarithmic terms are integrated up to infinity. 
In these last integrals, however, we assume that NJL interactions vanish for energies above the cutoff, therefore we 
set ct = A = for \p I > A . 

III. RESULTS 
A. Non-superconducting case 



As in most of the effective models of QCD interactions, the presence of a 2SC phase shows up at low temperatures 
and moderate chemical potentials. In order to study the effects of imposing the condition of color neutrality, let us 
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start by considering the NJL model without quark-quark interactions (H = 0), which means to set A = in the 
effective mean field thermodynamical potential in Eq. ([5]). Notice that, even in the absence of color superconductivity, 
one still may require a nonvanishing color chemical potential /ig owing to the presence of the quark interactions with 
the Polyakov loop. According to the discussion in the previous section, we take fi^ — 0, (j>s — 0, leading to a real- 
valued mean field thermodynamical potential. Thus the set of six relations ([6][7]) is reduced to only three nontrivial 
equations; two for the order parameters 

d f^MFA = Q (12) 



(da,d<t> 3 ) 

and one for the constraint of color charge neutrality 

<9^MFA 

dfi 8 



P8 = . (13) 



If we do not include quark-quark interactions, the model includes only three free parameters, namely the quark 
current mass mo, the coupling constant G and the cutoff A. For defmiteness, to perform the numerical analysis we 
will use a standard set of values for these parameters, taken from Refs. flol. [Tl|: 

mo = 5.5 MeV , G = 10.1 GeV~ 2 , A = 650 MeV . (14) 

We start by considering the case of T = 0. In this limit the Polyakov loop decouples, since the mean field value cf>3 
can be removed from the momentum integrals by a shift of the integration variable p^. One gets then two coupled 
equations, namely 



M-m M [ A , p 2 
= -7 / *7T 



G 7T 2 7 r E 



" 7^1 d "' 



2 Sg (£ + fi r ) + 2 Sg (£ - fi r ) + Sg{E + w ) + Sg [E - fi b ) 
Sg {E + fl r ) -Sg(E- fl r ) ~Sg{E + fl b ) + Sg(E- fl b ) 



(15) 
(16) 



For low values of the chemical potential fi, it is easy to see that both equations are trivially satisfied for a wide range 
of values of fig, while from Eq. (fTS"]) the value of M is fixed at M = Mo = 324.1 MeV. Thus, as expected, the system 
is in a phase in which the chiral symmetry is strongly broken. These results are valid as long as fi < Mo. The allowed 
values of fi% are those which satisfy \fi r \ < Mo and < Mq: 

Mp - fx < Ms < f (M + /i)/2 if /x<M /3 

2 V3 " \ Mo - fi if M /3 < fi < Mo 

Now if fi is increased above Mo the value of the effective mass M slightly decreases, while from Eq. (fTB|) one has 
fXg = 0. This holds up to a critical value fi c (T = 0) = 342 MeV. At this point one finds a first order transition into a 
normal quark matter (NQM) phase, in which the chiral symmetry is approximately restored. This behavior is shown 
by the T = curve in the upper panel of Fig. 1, where we plot the value of a = M — mo as function of the chemical 
potential fi. 

In order to determine the effect of the Polyakov loop one has to turn on the temperature, solving Eqs. (|12j|13j) . Let 
us first consider the region of low temperatures (T <C /i) and values of fi up to ~ Mo- In this region, while the value 
of the effective mass is kept at M ~ Mo, it is possible to show that the value of </>3 is basically determined by the 
minimization of the effective potential W($,T), which leads to 03 ~ 27r/3 (or, equivalently, $ ~ 0). In particular, it 
is interesting to consider the range < fi < Mo/3, where it is possible to perform some analytic calculations to solve 
Eq. (TT3")) . From this analysis, the value of fig is found to be given by the simple expression 

-^L = 2 it - Tln2 . (18) 
V3 

This result allows us to find the actual value of fig in the limit T = 0: one has ^g/v3 = 2/i, which leads to fi r = 3/i, 
fib = —3/i. Notice that the condition fi < Afo/3 is necessary to be in agreement with the constraints for fig, that we 
established above in the T = case. For low temperatures and -Mo/3 < fx < M , it is not easy to get a relation as 
simple as that in Eq. (fTS)) . However, it can be seen that in the T = limit Eq. (TT3^) implies fi r = Mo, which means 
/is/v3 = Mo — fi. These results are shown by the T — curve in the central panel of Fig. 1, where we quote the 
value of fig as function of the chemical potential. 



5 



In this way, it has been shown that for nonzero temperatures and baryon chemical potentials the interaction between 
the fcrmions and the Polyakov loop requires the presence of a nonzero color chemical potential ^ m order to keep 
color neutrality. This still holds in the T = limit if one demands continuity in the values of /is and fa. Hence the 
fact that in this limit the Polyakov loop decouples does not mean that the system is blind to the color symmetry 
breaking induced by fa. 

For finite values of the temperature, the critical values of the chemical potential /i c (T) can be found numerically, 
defining a first order phase transition line in the T — fi plane. This holds up to a given "end point" (/i£,Tg;), which 
is found to be located at about (333 MeV, 64 MeV) for our parameter set. For T > Tg the transition to the normal 
quark matter phase proceeds as a smooth crossover. The corresponding transition line can be established by looking 
at the peak in the chiral susceptibility, defined by 

1 d 2 ^ 

XS -2^?- (19) 

Our numerical results are shown in Figs. 1 and 2. In Fig. 2 we represent the T — fi phase diagram, denoting by solid 
and dashed lines the curves corresponding to first order and crossover chiral restoration transitions respectively. For 
comparison we also include the corresponding curves for the PNJL model without the imposition of color neutrality. 
It can be seen that the first order transition line extends up to about Tg = 100 MeV in this latter case. Notice that at 
fi = the phase transition temperature is found to be about 255 MeV, somewhat larger than the values obtained from 
lattice QCD |2|, [3j. This is a typical feature of local NJL models, and we do not adopt in our calculations the rescaling 
procedure 28] of the Polyakov-loop potential which leads to a flavor-dependent lowering of the critical temperature. 
It has been demonstrated that in nonlocal generalizations of the PNJL model [2!| the critical temperature at fi = 
is in agreement with the lattice result. In Fig. 1 we show the behavior of the effective mass a, the color chemical 
potential fis and the traced Polyakov loop $ as functions of the chemical potential, for some representative values 
of the temperature. The curves for T — 100 MeV and T = 200 MeV show clearly the crossover transition, while for 
T = 300 MeV the system is always in the normal quark matter phase (see Fig. 2). For relatively low values of fi, 
the central and lower panels of Fig. 1 show the correlation between the deconfincmcnt transition and the restoration 
of color symmetry: when the temperature is increased, the average mean field value of the traced Polyakov loop 
rises from (confinement) towards 1 (deconfincmcnt); this comes together with a decrease of the color chemical 
potential fig, which has been introduced in order to recover color neutrality. The effect is diluted in presence of a 
large baryon chemical potential. Finally, in the upper and lower panels we have also included with dashed lines the 
curves corresponding to the PNJL model without color neutrality (in the central panel, these correspond simply to 
fis — 0)- The results for the models with and without color neutrality are approximately coincident in the case of 
T = 300 MeV. 



B. Color-superconducting model 

Let us now take into account the quark-quark pairing interaction in Eq. ([T]), with the ansatz A5 = A7 = 0, A2 = A. 
According to the discussion in Sect. II, we have to solve the system of equations (fl"2"|) , together with a fourth equation 
<9J7mfa/<9A = 0. To proceed with the corresponding numerical calculations one has to fix the value of H, which 
cannot be directly obtained from T = fi = phenomenology. For definiteness we will consider values of the ratio H/G 
around 0.75, which is the value obtained from a Fierz rearrangement of the local quark-quark interaction in Eq. (jTJ) . 

Our main results are shown in Figs. [3] and [4] In Fig. [3] we show the arising T — fi phase diagrams. In order to see 
the qualitative dependence of the results on the value of H, we have considered the cases H/G = 0.7, 0.75 and 0.8. It 
is seen that now the phase diagrams include in general a large region of two-flavor superconducting (2SC) phases. For 
large values of the chemical potential and relatively low temperatures, the system is in a 2SC phase in which the chiral 
symmetry is approximately restored. Then, at T ~ 150 MeV one finds a second order phase transition into a normal 
quark matter (NQM) phase. In addition, for intermediate values of fi, the system gets into a 2SC phase in which 
the chiral symmetry is still strongly broken (%SB). Here the quarks acquire large dynamical masses, which are not 
substantially different from those obtained in the non-superconducting model. The size of this phase region is large 
for H/G > 0.75, and becomes reduced if the ratio is decreased. Once again, one finds a second order phase transition 
into the non-2SC phase, while the chiral restoration is driven by a first order phase transition for low temperatures, 
and by a smooth crossover for temperatures above a given end point. The shrinking of the xSB-2SC coexisting phase 
appears to be the main qualitative effect of the change in H/G within the considered range. 

In Fig. [4] we show the behavior of the mean field diquark condensate A (left panel) and the traced Polyakov loop $ 
(right panel) as functions of the chemical potential for T = 0, 50 MeV, 100 MeV and 150 MeV. The curves correspond 
to the case H/G = 0.75. Both the first order and the crossover chiral restoration transitions can be noticed through 
their effects in the diquark condensate, while the size of the diquark gap is found to be of the order of a hundred 
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FIG. 1: Behavior of the effective mass a, the color chemical potential and the traced Polyakov loop $ as functions of the 
chemical potential, for various values of the temperature, in the non-superconducting two-flavor PNJL. Solid and dashed lines 
correspond to the results with and without the imposition of color neutrality, respectively. Notice that in the upper and lower 
plots the solid and dashed curves corresponding to T — 300 MeV are almost coincident. 



MeV for the bulk of the temperature range (it decreases rapidly to zero near the phase border, at T ~ 150 MeV). 
Concerning the traced Polyakov loop, it is seen that the behavior is similar to that found in the case of the non- 
superconducting model, in the sense that $ < 0.15 for temperatures below 100 MeV and chemical potentials up to 
500 MeV (see lower panel of Fig. 1). Here the rise with \x is even less pronounced than in the non-superconducting 
case, since one finds a plateau in the xSB-2SC coexistence region (see curves corresponding to T = 50 and 100 MeV). 
For comparison, the behavior of 4> in the non-superconducting model for T = 100 MeV has also been plotted in the 
figure. For T > 150 MeV there is no 2SC phase, thus there is no difference between the behavior of $ in both models. 
Finally, the behavior of the color chemical potential /is i n both superconducting and non-superconducting models is 
also found to be qualitatively similar (see central panel of Fig.[T]). 

It has been noted earlier by Fukushima [2(| that the phase diagram of the PNJL model exhibits at densities \i > \i c 
and low temperatures a so-called quarkyonic phase characterized by the coexistence of confinement ($ « 1) and 
chiral symmetry restoration (see also [22j). The possibility of such a phase was suggested in Ref. [30( for large N c 
QCD where also the term was coined. In a different context it has been discussed in Ref. [3l|. It is at present an 
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FIG. 2: Phase diagram for the non-superconducting two-flavor PNJL model with and without the imposition of color neutrality. 
Solid and dashed lines correspond to first order phase transition and smooth crossover, respectively. The fat dots denote the 
position of the end points. 
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FIG. 3: T — fi phase diagrams for the neutral superconducting two-flavor PNJL model, for various values of the ratio H/G. 
Solid, dashed and dotted lines denote first order, second order and crossover phase transitions, respectively. The fat dot 
indicates the end point. 




FIG. 4: Behavior of the color superconducting condensate A (left) and the traced Polyakov loop $ (right) as functions of the 
chemical potential, for various values of the temperature and a coupling ratio H/G — 0.75. For comparison, the dashed-dotted 
line in the right panel shows the curve corresponding to T — 100 MeV in the non-superconducting model (H — 0). 
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open question whether a quarkyonic phase shall exist in the real world with N c — 3. As a result of the present study, 
we find that the requirement of color neutrality together with the occurrence of diquark condensation are consistent 
with the existence of a quarkyonic phase in the PNJL model. Moreover, as it is shown in the right panel of Fig. |4j 
the presence of diquark condensation leads to a reduction of the expectation value of the traced Polyakov loop <E> and 
thus enforce the confining aspects of the model. 

It is important to point out that the presence of a nonvanishing Polyakov loop variable fa breaks the color SU(3) 
symmetry down to SU(2), so that the rotational invariance for the choice of the orientation of the 2SC gap under color 
neutrality constraints 32j might be lost. As a consequence, the above ansatz for A Q may not correspond to the true 
minimum of the thermodynamical potential which instead might point to a different direction in the A a space [33]] . 
While accounting for such a disorientation of the 2SC condensate may lead to quantitative corrections, we believe 
that the qualitative structure of the phase diagram for the PNJL model under color neutrality constraints shown in 
the left panel of Fig. [3] remains robust. The 2SC phase might coexist with chiral symmetry breaking and confinement 
in the temperature range between about 50 MeV and 100 MeV at densities below the onset of chiral restoration with 
possible observable consequences in heavy-ion collision experiments, e.g., due to diquark- antidiquark annihilation into 
lepton pairs (34[ . Note that with the coexistence of %SB and 2SC phases the condition for Bose-Einstein condensation 
(BEC) of diquark bound states is fulfilled and the chiral restoration transition in the above temperature window may 
therefore be characterized as a BEC-BCS crossover transition where the diquark bound states undergo a Mott effect 
to unbound, but resonant diquark Cooper pairs [35L l36j 



IV. CONCLUSIONS 



In this paper we have explored the consequences of imposing the color neutrality constraint within the framework of 
the standard mean field treatment of the PNJL model at finite chemical potential. We have found that the fulfillment 
of such a constraint leads to some changes in the phase diagram of the model, as e.g. a shrinking of the xSB domain 
in the T — /i plane and a shift of the end point to lower temperatures. The most noticeable effect turns out to be the 
presence of a coexistence region of %SB and 2SC phases, which may have consequences for possible applications of 
the model in heavy-ion collision experiments. 

The effects can be understood due to the nonvanishing color chemical potential which is introduced to com- 
pensate the color charge ps induced by the color background field fa. For a given baryochemical potential /i, fi$ > 
reduces the chemical potential of blue quarks, while simultaneously increasing that of red and green ones. The net 
effect is a shift of the x^B phase border to lower values than without the color neutrality constraint, due to the 
nonlinear dependence of the quark mass gap on the chemical potential close to its critical value for a given, not too 
low temperature. It is seen that this effect is relatively slight, taking into account the theoretical uncertainty that 
can be in general expected from an effective quark model. In addition, the increase of chemical potentials of red and 
green quarks, which pair in the 2SC channel, entails an effective lowering of the critical value of fi = /i c for the onset 
of the 2SC phase, again for not too low temperatures. One finds in this way a xSB-2SC coexisting phase, whose size 
depends on the value of the parameter H that drives the quark-quark interaction. We notice that the 2SC border has 
been obtained after choosing a given ansatz for the mean field values A a . The actual orientation of the 2SC gap in 
color space under the presence of the Polyakov loop is a subject that deserves further study. Moreover, it is worth to 
point out that at finite densities there is no fundamental need to be restricted to Lorentz invariant condensates, thus 
more general mean field configurations could be considered. 

Signatures of the coexisting %SB and 2SC phases might become apparent in not too energetic heavy-ion collision 
experiments (about 3-5 GeV/nucleon), as they are possible at the future facilities FAIR-CBM (Darmstadt), NICA- 
MPD (Dubna) and J-PARC (Japan). Such a mixed phase can give rise, e.g., to an enhancement of low-mass dilepton 
production due to diquark- antidiquark annihilation. It is a realization of the BEC-BCS crossover in quark matter. 
While such a transition has previously been discussed in the NJL model by artificially increasing the diquark coupling 
to large values [H, H?],[38j], which can result in conflicts with the nuclear matter ground state, in the present PNJL 
model the xSB-2SC coexisting phase has been obtained without any artificial tuning of parameters. It just occurs as 
a consequence of color neutrality in a domain of temperatures and chemical potentials around the critical end point 
of the first order chiral phase transitions for the standard Fierz- value of the diquark coupling, H/G = 0.75, without 
affecting the nuclear matter ground state at T = 0. 

Finally, it has been shown that the coexistence of confinement and chiral symmetry restoration characterizing 
a quarkyonic phase in the PNJL model at finite densities and low temperatures gets even reinforced due to color 
superconductivity. The occurrence of diquark condensation reduces the expectation value of the traced Polyakov loop 
and thus acts towards strengthening the confinement aspect of the model. 
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